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Context. The largest uncertainty for cosmological studies using clusters of galaxies is introduced by our limited knowledge of the 
statistics of galaxy cluster structure, and of the scaling relations between observables and cluster mass. A large effort is therefore 
undertaken to compile global galaxy cluster properties in particular obtained through X-ray observations and to study their scaling 
relations. However, the scaling schemes used in the literature differ. 

Aims. The present paper aims to clarify this situation by providing a thorough review of the scaling laws within the standard model 
of large-scale structure growth and to discus various steps of practical approximations. 

Methods. We derive the scaling laws for X-ray observables and cluster mass within the pure gravitational structure growth scenario. 
Using N-body simulations we test the recent formation approximation used in earlier analytic approaches which involves a redshift 
dependent overdensity parameter. We find this approximation less precise than the use of a fiducial radius based on a fixed overdensity 
with respect to critical density. 

Results. Inspired by the comparison of the predicted scaling relations with observations we propose a first order modification of the 
scaling scheme to include the observed effects of hydrodynamics in structure formation. This modification involves a cluster mass 
dependent gas mass fraction. We also discuss the observational results of the reshift evolution of the most important scaling relations 
and find that also a redshift dependence of the gas mass to total mass relation has to be invoked within our modification scheme. 
Conclusions. We find that the current observational data are within their uncertainties consistent with the proposed modified scaling 
laws. 
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1. Introduction 

Galaxy clusters form from overdense regions in the large-scale 
matter distribution, which have small amplitudes at early epochs 
and only collapsed to objects very recently. In the standard cos- 
mological model the large-scale matter distribution is described 
by a random Gaussian field characterized by a power spectrum 
with a smoothly changing power law index over relevant length 
scales. This implies that the structure evolution will feature a 
large degree of self-similarity in scale and time (e.g. Peebles 
1980). Consequently galaxy clusters, which form an integral part 
of this large-scale structure, also show an imprint of this general 
self-similarity. This connection between the framework of the 
evolution of the gravitating matter on large scales and galaxy 
cluster formation and their observed structure was realized in 
early studies, e.g. by Gunn & Gott (1972), Fillmore & Goldreich 
(1984), Bertschinger (1985), Hoffman & Shaham (1985) and in 
simulations e.g. Frenk et al. (1985), Zureck, Quinn & Solomon 
1988, Efstathiou et al. 1988, West, Dekel & Oemler 1987). It 
resulted in a comprehensive description of the structure of dark 
matter halos, of which galaxy clusters are the most massive rep- 
resentatives, in a series of papers by Navarro, Frenk & White 
(1995, 1996, 1997) and follow-up literature. In this picture of 
purely dark matter structure growth, dark matter halos form a 
nearly self-similar, two-parameter family, with the two param- 
eters being mass and a concentration or time-of-formation pa- 
rameter. This structural model of clusters describes an average 
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behavior of the cluster population, where the different statistical 
realizations of mass distributions in the protoclusters produce a 
significant scatter in the observed structural parameters around 
this mean. Deviations from the equilibrium state after merger 
events further contribute to this scatter. 

The addition of baryons to this model leads to a modification 
of this picture, which for clusters can be seen as a perturbation of 
the dark matter structure evolution. In this sense galaxy clusters 
mark the very interesting transition region, where a first order 
description involving only the large-scale structure gravitational 
physics provides a very effective guideline and the more compli- 
cated hydrodynamics, including radiative cooling and feedback 
from star formation as well as AGN activity, constitutes a pertur- 
bative refinement. At smaller scales, for galaxies the gaseous as- 
trophysics acting on small scales becomes dominant for the ap- 
pearance of the visible objects and the observed evolution of the 
large-scale structure on galaxy scales becomes very non-linear. 
It is therefore on galaxy cluster scales where we can still very 
successfully apply analytical descriptions as a useful guideline 
for the understanding of structure evolution. The paper builds on 
this property of galaxy clusters. 

X-ray observations are currently providing the most detailed 
account of galaxy cluster structure and are consequently used 
extensively to test the predictions of the large-scale structure 
growth models. However, they do not directly provide a picture 
of the dark matter halo distribution, but the distribution of the 
hot intracluster medium (ICM) that fills the entire cluster vol- 
ume and radiates in X-rays. Therefore the X-ray appearance of 
clusters includes aspects of the hydrodynamics how the gas re- 
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acts to the dark matter density distribution and how the ICM 
evolves in its thermodynamic properties (e.g. Voit 2005). One 
can use it as a tracer of the dark matter distribution, for example 
through the assumption that it is located in the dark matter po- 
tentials in hydrostatic equilibrium. We can thus expect that the 
X-ray appearance of galaxy clusters is featuring some hydrody- 
namic modification compared to the more readily described dark 
matter distribution. 

Because the galaxy cluster formation is so tightly connected 
to the large-scale structure evolution and the fact that there is a 
well described self-similar structure statistics in the purely gravi- 
tational cluster formation model (Navarro et al. 1997) we can ex- 
pect that there are simple analytically derivable scaling relations 
for the global X-ray observables as a function of cluster mass. 
These relation have been studied already early by e.g. Kaiser 
(1986), and Evrard & Henry (1991) and this theoretical work 
has been supported by simulations (e.g. Bryan & Norman 1998, 
Borgani 2004, Kravtsov et al. 2006, Evrard et al. 2008, Stanek 
et al. 2010, Short et al. 2010, Borgani & Kravtsov 2010). 

With the event of detailed observational studies of cluster 
structure in X-rays by means of the advanced X-ray observato- 
ries Chandra and XMM-Newton, large sets of observational data 
on cluster structure and scaling relations have become available 
now and the detailed testing of the theoretical predictions for the 
scaling laws is in full swing (e.g. Markevitch et al. 1998, Arnaud 
& Evrard 1999, Mohr & Evrard 1997, Finoguenov et al. 2001, 
Ikebe et al. 2002, Reiprich & Bohringer 2002, Ponman et al. 
2003, Ettori et al. 2004, Vikhlinin et al. 2005, Pointecouteau et 
al. 2005, Arnaud et al. 2005, Pratt et al. 2006, Kotov & Vikhlinin 
2006, Zhang et al. 2006, Maughan et al. 2006, Maughan 2007, 
Arnaud et al. 2007, Pratt et al. 2009, Mantz et al. 2010, Arnaud 
et al. 2010, Sun et al. 2011, Reichert et al. 2011). An investi- 
gation of the relevant literature shows, however, that a number 
different methods are used for the scaling of the data at different 
redshifts. The aim of this paper is therefore to critically review 
these methods and to determine the best approach based on com- 
parison with simulations and observations. 

In the above mentioned literature mostly analytical formula- 
tions of the scaling relations have been used, based on general 
considerations of structure formation. In order to provide the 
ground for higher precision in the analysis of the evolution of 
cluster structure, the simplifications made in the analytical mod- 
els should be replaced by tests and calibrations with N-body sim- 
ulations. This situation can be compared to that of the theoretical 
prediction of the dark matter halo (galaxy cluster) mass function, 
where the analytical model by Press and Schechter (1974) has 
paved the way for the general formulation of the solution for the 
mass function, but the actual formulae now applied, are the re- 
sult of careful calibration with N-body simulations (e.g. Jenkins 
et al. 2001, Evrard et al. 2002, Warren et al. 2006, Tinker et al. 
2008). Here we adopt a similar approach. We first present the 
theoretical framework for the description of the evolution of the 
scaling relations in the classical form based on the assumption of 
the recent formation approximation and compare it to an alter- 
native scheme used in the literature. We then resort to the results 
of N-body simulations to test the predictions of the analytical 
approaches and discuss which of the presently used methods in 
the literature is best. 

In the second part of the paper we compare the theoretical 
scaling relations to observations, inspect the deviations of the 
observed scaling relations from the predictions based on dark 
matter structure evolution (often called "gravitational scaling re- 
lations"), and discuss these deviations in the context of the in- 
fluence of hydrodynamical processes. We then explore a simply 



empirical modification scheme of the scaling relation using a 
mass dependent depletion factor for the ICM gas to account for 
the hydrodynamical scaling effects and compare the so obtained 
set of scaling relations to observations. In a last step we consider 
how this modification should depend on redshift to be consistent 
with the observational data. 

The paper is structured as follows. In section 2 we derive the 
"gravitational scaling relations" based on dark matter structure 
evolution with the assumption that the baryonic matter follows 
the dark matter. In section 3 we investigate the dependence of 
these scaling relations on the redshift dependent overdensity pa- 
rameter in a ACDM cosmology and test this model against the 
method using a fixed overdensity parameter in section 4 using 
numerical simulations. In section 5 we discuss the redshift de- 
pendence of the overdensity parameter, which defines the proper 
fiducial radius of the clusters, in the context of numerical stud- 
ies of the redshift dependence of the concentration parameter of 
galaxy clusters. Section 6 then starts the second part of the pa- 
per where we discuss the modification of the scaling relations to 
include hydrodynamical effects. After a comparison of the de- 
scription of the evolution of the scaling relations in terms of the 
parameters E(z) and (1+z) in section 7 and a brief comparison 
with some recent simulations in section 8, we provide a compre- 
hensive comparison of scaling relation results in the literature 
with the model predictions in section 9. Finally section 10 con- 
tains a discussion and conclusions. 



2. Analytic gravitational self-similar model 

In describing the self-similar structure of galaxy clusters in the 
purely gravitational picture, we will first consider the cluster for- 
mation by collapse in an Einstein-deSitter model (EdS), that is a 
Universe with critical density and zero cosmological parameter. 
When we describe the self-similar evolution for other cosmolo- 
gies in the second step, this model is used as reference model. 
For illustration, as sketched in Fig. 1, we use a "top-hat overden- 
sity" for the initial conditions characterized by a homogeneous 
overdensity within a sphere. Under realistic conditions proto- 
clusters will have a wide range of morphologies, but we can rea- 
sonably assume, that the morphology distribution with respect 
to the top-hat model is similar for different masses or formation 
times. The statistical realizations are strictly self-similar only as 
far as the power spectrum of the density fluctuations, P(k), is 
described by a power law. 

The theoretical background to the proper scaling has been 
worked out some time ago, with a seminal paper being provided 
for example by Kitayama & Suto (1996). The aspects of this 
model, which are crucial for our discussion, are illustrated in 
Fig. 1, and will be first interpreted in the frame of an EdS cos- 
mological model. The left hand side of Fig. 1 shows a cluster 
which virializes at redshift zero, where the point of virialization 
is in general defined by the time when an ideal homogeneous 
sphere would have collapsed to a point. A good time of refer- 
ence in the collapse process is the epoch when the overdensity 
stops expanding and turns around to collapse, which happens 
for this cluster at z=0.78 when the mean matter density inside 
the protocluster is about 6.5 times the background density. This 
is true for all clusters with different mass which virialize at z=0. 
All these clusters finish their formation at the same time - twice 
the turn-around time - and in the same way from the spherical 
overdensity with same amplitude to a self-similar structure with 
same density shape and amplitude just with a different radial 
scaling. Therefore, for clusters formed at the same epoch we can 
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Fig. 1. Schematical illustration how the galaxy cluster structure 
scaling depends on the formation epoch of the cluster and on the 
critical density of the Universe in an Einstein-de Sitter (densi- 
ties have blue color) and in a concordance cosmological model 
(densities have green color). The left panel provides a sketch of 
cluster collapse for a cluster with a formation redshift of z — 
and the right panel shows a cluster forming at z — 0.5. For further 
explanations see the text. 



defined a fiducial outer radius by a mean cluster density thresh- 
old in units of the mean or critical density of the EdS universe. 

Inspecting now the evolution of a cluster that forms at higher 
redshift, like the cluster sketched on the right hand side of 
Fig. 1, we find a similar picture. The protocluster has a higher 
mean density when it forms compared to younger clusters, but 
the ratio of mean protocluster density to the background den- 
sity at turn-around is the same factor of ~ 6.5. From the turn- 
around epoch, cluster collapse will be self-similar to the forma- 
tion at earlier (later) epochs, just scaled to higher (lower) density. 
Also the background density of the universe will evolve self- 
similarly during cluster collapse, since in the EdS universe we 
have p,„ = p cr i, oc r 2 . Therefore, if we define a fiducial outer ra- 
dius by the radius at which the mean density of the cluster has the 
same overdensity ratio, A, to the mean density of the universe, 
we can compare self-similar radii for clusters at different epochs 
and different sizes. For the EdS a popular choice for the fiducial 
radius ,r&, is for example the virial radius (with A ~ 18 n 2 ~ 180; 
e.g. Peebles 1980). 

A self-similar fiducial radius of the clusters is thus defined 

by: 



? - 

' A — 



Anp crit h 



M c luster(r < r A ) 



(1) 



The picture becomes more complicated when we change to a 
low density universe (possibly also including a A term). For the 



overdense regions developing into clusters there is no change. 
According to the Birkhoff theorem, where a local region of the 
universe evolves like a universe with these local density and ex- 
pansion parameters irrespective of the embedding cosmology, 
the cluster evolution does not care about the background uni- 
verse 1 . Therefore we can keep the knowledge we have gained 
in the EdS reference frame, we just have to introduce another 
conversion which links the critical density of a general model 
universe with the density of a coevolving EdS universe. To make 
this more transparent, we illustrate this point for the case of the 
Concordance Cosmological Model (CCM, with Q„, = 0.3 and 
Q.A = 0.7) which seems to provide a close approximation to the 
structure of the real Universe (e.g. Spergel et al. 2007, Komatsu 
etal. 2011). 

This is again illustrated in Fig. 1 where a lookback time of 
13 Gyr (corresponding to a redshift of z=10 in the CCM) was 
used as a good approximation to set an EdS and CCM model to 
the same density initial conditions and look at their coevolution 
in time which is shown in Fig. 2. Note that for a given time the 
redshifts of the two models will differ. We see that approaching 
present time the CCM model has a lower density than the EdS 
model, but a higher critical density, which is due to the accelerat- 
ing expansion. Defining the ratios of the critical densities in Fig. 
2 as retime) = t](z C cm) = p cr it(CCM)/p m (EdS), where z C cm is 
the redshift in the CCM (with an index that we will drop in the 
following), we can adapt the overdensity parameter to the new 
situation by defining 



A(Xm(z) = &EdS h(z) 
and 

^ccm(z) = A C cm(z = 0) * r](z = 0)/tj(z) 



(2) 



(3) 



Thus, as for example illustrated in Fig. 1, a Asds - 200 in 
EdS corresponds to a Aqcm ~ 101 for z — in the CCM model 
since jj(z = 0) ~ 2. It corresponds to A ~ 138 at z = 0.5 in the 
CCM model, however, since with increasing redshift the differ- 
ence between the two models shrinks. In the following we drop 
the index of A and use it only for the CCM. This formalism has 
been worked out in detail for the general case by Kitayama & 
Suto (1996), Eke et al. (1998) and useful approximate formulae 
for A(z) for a range of cosmological models including non-flat 
universes can e.g. be found in Pierpaoli et al. (2001). 

For the dependence of the fiducial radius on mass we have 
a simple geometrical scaling, while for the scaling with time we 
have a proportionality of the mean density of the cluster to the 
critical density of the universe, po, taken either at the time of 
turn-around or at collapse (for the EdS model). The mean density 
as a function of redshift is: 



Pcritical(z) 
Pcriticaliz = 0) 



H(z) 2 
Hi 



= E(zf 



(4) 



The evolution of the radius at fixed overdensity (e.g. A = 
200) is thus given by 



I M 20 q 
\ Po 



1/3 



(5) 



As mentioned above, the radius at fixed overdensity de- 
scribes a self-similar region in clusters at different epochs only 



1 The cosmological constant introduces a slight change in the col- 
lapse evolution of the order of a percent, which is neglected here 
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Fig. 2. Evolution with lookback time of the mean density, p,„ and 
critical density p c ,„ in the Concordance Cosmological Model in 
comparison with the matter density of a coevolving Einstein-de 
Sitter reference model, po. The curves are from top to bottom: 
Pent, Po, and p,„. The two models are compared at the same times 
after the starting point of the coevolving calculations at 13 Gyrs. 



in the EdS scenario. Thus if we keep true to the described mod- 
eling, and if we want to compare like with like, we have to use 
the radii, r&( z ) for any comparison. Thus Eq. 5 becomes 

r m oc M^£(z)- 2/3 A(z)- 1/3 (6) 

We ought to note here, that in this approach it has been as- 
sumed, that the clusters we observe have just collapsed. This is 
true in a very broad sense only, since clusters are always ac- 
creting material and we always see them in a stage where they 
have just completed some late accretion. This scenario is termed 
"recent formation approximation" in the literature. A more crit- 
ical inspection of cluster formation e.g. in N-body simulations 
shows, however, that the cluster structure depends in more detail 
on the whole recent accretion history. We should therefore be es- 
pecially concerned about this approximation within the frame of 
the presently prefered ACDM cosmology. While in an Einstein- 
de Sitter cosmological model structure growth is continuing into 
the future, it starts to cease in a universe with a low matter den- 
sity as soon as the matter density drops below the critical density 
value. Thus clusters have accreted matter more slowly in the re- 
cent past than at higher redshift and will practically stop to grow 
in the distant future (e.g. Busha et al. 2007). Thus even though 
we recognize the beauty and logic of the above approach, we 
have to critically test how much deviations are introduced by 
these approximations by finally comparing to N-body simula- 
tions. 

After the radius - mass relation given in Eq. 4, the next basic 
equation is the one linking the X-ray gas temperature with the 
cluster mass. The heat of the ICM comes from the conversion of 
potential energy during the formation of the cluster. We there- 
fore expect the temperature to be proportional to the depth of the 
gravitational potential and we thus find: 

T X cc Co cc M^Etf* A(z) 113 (7) 

In a similar way other essential scaling relations for impor- 
tant X-ray properties can be constructed (see also Kaiser 1986). 



In Table 1 we list a set of scaling relations involving the cluster 
mass or the temperature as scaling parameter (the latter being a 
prime observable parameter) for X-ray determined properties as 
X-ray luminosity, Lx 2 , gas mass, M gas , ICM entropy, K, den- 
sity, p, Tx-parameter, with Yx = M gas x T x , and pressure, P. We 
have termed these relations "gravitational", since they only in- 
clude the physics of structure evolution of the dark matter which 
only interacts gravitationally. The baryonic matter is assumed 
here to follow the dark matter and hydrodynamical effects have 
been neglected. Consequently we will term the relations taking 
these effects into account as "hydrodynamical relations", which 
better describe the observations (e.g. Voit 2005). 

In these relations the redshift scaling factor E{z) 2 x A(z) ap- 
pears always in the same combination just with different pow- 
ers, and thus one can introduce the abbreviation, F(z) = E(z) x 
A(z) 1 ^ 2 . Note that all the integral observables or parameters ap- 
pearing in Table 1, like Lx, M tot , M gas , Yx have to be integrated 
out to /?a(z)> if a comparison between clusters at different red- 
shifts are made. Thus, M in the table should actually be written 
as Ma(j); we have chosen the simplified version of the formulae 
in the Table for easier reading. 

3. Dependence on the overdensity parameter 

The appearance of the parameter A(z) in the above equations is a 
nuisance, in particular as it can only be calculated by numerical 
integration or from approximate formulae given in the literature, 
e.g. Pierpaoli et al. (2001). Thus some effort has been done to 
check with simulations if this parameter is really necessary, and 
if the scaling relations can simply be derived from simulations 
without a redshift dependent overdensity parameter. The work 
by Evrard et al. (2002) and (2008) is for example performed in 
this spirit. 

On the other hand we have seen in the previous chapter, that 
this parameter has nothing to do with the cluster formation. It 
comes solely from the break in the self-similar evolution in the 
background cosmology in going from an EdS to a CCM sce- 
nario. Therefore the overdensity parameter should not be easily 
abandoned without further checking. 

In the following we will investigate if the unwanted param- 
eter can be eliminated if we make assumptions on the structural 
parameters of the clusters. One of the most fundamental rela- 
tions is that between temperature and mass. Making the very 
simple assumption that the mass profile is given by an isothermal 
sphere, where M(r) cc r, we find with some arithmetics involving 
Eq. 1 that M A oc A~ 1/2 . 

Eq. 7 can be rewritten as 

M m cc Tf~ E(z)- 1 A(zr 1/2 (8) 
with the above suggested mass dependence on A(z) we find 

that 

M A(/ , f) (z) = M A( ^^j^J ; <x7-f £(z)-'A(/«r l/2 (9) 

where A(fix) is the fixed value for the A parameter inde- 
pendent of redshift. We note that for the case of the isothermal 

2 For the scaling of the X-ray luminosity in a specific band, L x , we 
assume that the energy band is chosen such, that the X-ray emissivity 
is independent of the ICM temperature. This is for example almost ful- 
filled for the 0.5 to 2 keV energy band used in most imaging analysis 
of galaxy clusters, where the change in emissivity in the temperature 
range from 2 to 10 keV is less than 6% for given emission measure. 
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property 


proportionality 


scaling for var.A 


r - dependence' 1 


scaling for fixed A 


deviation 


radius,,^ 


oc M 1/3 F- 213 


oc T 1 ' 2 F- 1 


cc r 


T U2 £ -l 




density profile" 


p cc M if" 3 


oc F(z) 2 


oc r 


oc E(z) 2 


- 


luminosity'' 


L x oc p 2 R 3 


cc r 15 f 


oc r 1/2 


oc T| 5 £ A(z) 3/4 


39% 


bolom.lum. 


L bol oc p 2 R? T 1 ' 2 


cc T 2 F 


oc r 1 ' 2 


cc T 2 E A(z) 3/4 


39% 


entropy, K 


oc T/p 2 ' 3 


cc T F~ 4/3 


o<- T 


oc r x £- 4/3 A( z )- 1/6 


7.5% 


gas mass 


ccpR 3 


cc T 3n F _I 


oc /■ 


oc r 3/2 £-' 


- 


Y x 


cc TpR 3 


cc T 5 ' 2 F- [ 


oc r 


oc T^E- 1 


- 


pressure 


oc T p 


cc r f 2 


oc r- 7/3 


oc r x £ 2 A(z)- 1/6 


7.5% 


surf.bright/' 


Sx K p 2 R 


cc T" 2 F 3 


oc r~ 3 


oc T\l 2 E 3 


- 


temperature 


T x 


oc M 2 ' 3 F 2 ' 3 


(M(r) oc r) 


oc M 2/3 £ 2 ' 3 




luminosity* 


L x 


ocMF 2 




oc M+E 2 A(z) 3/4 


39% 


bolom.lum. 


Ebol 


oc M 4 ' 3 F lj3 




cc M 4 / 3 E 1 ' 3 A(z) 3 ' 4 


39% 


entropy 


K 


oc M 2 ' 3 F- 2 ' 3 




oc M 2/3 £- 2 ' 3 A(z)- 1 ' 6 


7.5% 


gas mass 




oc M 




oc M+ 




Y x 


Yx 


oc M 5 ' 3 F 2 ' 3 




cc £ 2 ' 3 




pressure 


P 


oc M 2 ' 3 F g/3 




cc M 2 J 3 E*' 3 A(zr" 6 


7.5% 


surf.bright." 


Sx 


oc M [/3 F [0/3 




cc M [ J 3 E 10 ' 3 





Columns 1 and 2 give the property and its definition, column 3 provides the scaling relation where the parameter F(z) can be read in two 
ways: (i) for the recent formation approximation approach all quantities which involve a radial integration have to be taken inside r & and 
F(z) = E(z) x A(z) 1 ' 2 (note that the property of column 1 is for this case assumed to be taken at the radius of the overdensity A(z) ), 
alternatively it can be used for the fixed overdensity model, then all the radially dependent properties have to be taken at r A and F = E(z). 
Column 4 gives the assumed radial dependence of the property of column 1 and column 5 the scaling relation for the case that the property 
of column 1 is taken at a fixed overdensity radius. The temperature T x is in all cases defined as a mean measured temperature, determined in 
practical terms usually as a mean temperature in the region between two defined overdensity radii, except for entropy, K and pressure, where it 
means the temperature at the fiducial radius. In general luminosity, gas mass and Y x are integral parameters, while density, entropy, pressure, 
and surface brightness are defined as local parameters at the fiducial radius. M* is the mass at fixed overdensity. Column 6 lists the deviation 
of the scaling relation if the A(z) term is neglected for a comparison of clusters at z = and z = 1 . 

a the columns 4 to 6 refer to the radial dependence of the profile in the radial range r ~ rsoo to nooo 
Here we assume p oc r~ 2 , a y0-model for the gas with yS = 2/3 and K oc r which implies T oc r~ [,i 

b assuming a temperature independent emissivity for the luminosity in a restricted soft X-ray band (see text) 



sphere mass profile the redshift dependent overdensity parameter 
is eliminated from the equation of the mass-temperature relation. 
In the more general case of M(r) oc r y we find 

/200\^ 

r A = r 20 o — (10) 
and 

M A ocA^7 (11) 

The most widely used model for cluster mass profiles estab- 
lished by simulations and best confirmed by observations is the 
NFW model (Navarro et al. 1995, 1997). 3 In Fig. 3 we show 
the logarithmic slope of this mass profile as a function of ra- 
dius and the logarithmic slope for the function M(A(z)) for a 
typical value of c = 5. For the range of interest for overdensi- 
ties of 100 to 2500 the slope parameter for the latter function is 
in the range 0.3 to 0.6 not far off from the case of the isother- 
mal model and therefore we can expect an approximate elimina- 
tion of the A(z) parameter. Most of the observational results have 
been obtained for radii corresponding to overdensities of 500 for 

3 While recent N-body simulations show deviations from a NFW pro- 
file and promote improvements in form of e.g. an Einasto profile (e.g. 
Navarro et al. 2004, Gao et al. 2008), the NFW model is sufficiently 
accurate for our purpose. 



the CCM model which corresponds roughly to an overdensity of 
1000 for EdS at redshift zero and a smaller value for higher red- 
shifts. The isothermal sphere approximation best applies in the 
interval 500 to 1000 where the slope parameter for the M(A(z)) 
function is in the range 0.4 to 0.5. The maximum deviation in- 
troduced in the comparison between z — and z — 1 clusters at 
these overdensities is of the order of 4% and thus far smaller than 
the uncertainties in all observed relations and also smaller than 
the errors in relations derived from simulations. Thus, for the 
current precision of the results we can easily neglect the A(z) de- 
pendence in the mass-temperature relation. This will, however, 
not be true for any relation e.g. like those listed in Table 1. In 
this Table we apply the conversion to A(fix) formulation in a 
way analogous as done for the mass profile above, with results 
shown in column 5. For parameters like luminosity, entropy and 
pressure the A dependence does not cancel for the case of the 
recent formation approximation. 

The X-ray determined temperature and the velocity disper- 
sion of the galaxies in optical observations have generally been 
used as two of the most reliable proxies for the estimate of clus- 
ter masses 4 (e.g. Arnaud & Evrard 1999, Carlberg et al. 1996, 
Biviano et al. 2006). Therefore it is interesting that it is ex- 
actly for this relation, that the A(z) dependence can be neglected. 



4 In the recent literature notably in Kravtsov et al. 2006 the parameter, 
Y x , is often promoted as the best mass proxy. 
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Fig. 3. upper panel: logarithmic slope of the NFW mass profile 
as a function of scaled radius. Lower panel: logarithmic slope 
of the NFW mass profile a function of the mean density of the 
cluster over the critical density. 



Indeed in the observational work by Kotov & Vikhlinin (2006) 
for example the evolution of the mass-temperature relation is ex- 
plained without the need of A(z), and similarly in the simulations 
by Evrard et al. (2008) they detect a perfect relation for the one- 
dimensional dark matter velocity dispersion and cluster mass of 
the form 



M 2 „„occ^ 75±0 - 023 EizT 1 



(12) 



From these two examples we cannot conclude, however, that 
A(z) is an unnecessary parameter in general. This will be illus- 
trated through the scaling of density and radius with redshift. 
According to the relations given in Table 1 and from the illustra- 
tion in Fig. 1, the density scales as p(z) oc E(z) 2 A(z). To see how 
this works we need a reference density within the density profile 
of the cluster. In a cored density profile, we can use the central 
density, po- For the NFW density profile, described by 



Pcrit 5c 



(r/r s )(l+r/r s y 



(13) 



where p„„ is a reference density that is obviously propor- 
tional to A(z), we can for example use p s = p(r = r s ) as ref- 
erence density for comparison. This density is characterized by 



the break of the slope in the density profile. If we compare this 
density point in different clusters it should have the same over- 
density in the EdS cosmology and in the CCM it should scale 
according to the above relation. The radius of this point should 
then scale as given by Eq. 6 and expressed as a function of tem- 
perature like 



r A(z) (x T m E~ l A(z) 



(14) 



in which the A(z) dependence can be eliminated according to 
the relations used in Eq. 9 to yield: 



r Mfix) 



oc T 



1/2 



(15) 



This then implies that in comparing density profiles we will 
not observe a dependence on the A(z) parameter in the radial 
scaling around r^o to riooo, but we should observe the influence 
of this parameter in the amplitude scaling of the density profiles 
at smaller radii according to Table 1 : 



p oc E(zf A(z) oc Pcrit A(z) 



(16) 



From z=0 to z=0.5 (z=l) the parameter A(z) changes by a 
factor of 1.36 (1.55). We will test this relation in the next section. 

For any observable, Obs&, the influence of the A(z) parame- 
ter can be investigated in the following way (using a logarithmic 
Taylor expansion): 



Obs A ( z) = Obs A(fLx) 



dlnObs(r) dlnr 

A(z) dlnr dlnA 



A(fix) 



(17) 



where the last factor in the exponent is again approximately 
-0.5. For example for the total bolometric luminosity inside r& 
we find for a /?-model surface brightness profile with (3 ~ 2/3 a 
behavior of Ltoi(r) oc r° 5 and consequently 



L bo ,(A(z)) = L bo ,(A(fix) 



A(z) 
A(fix) 



-1/4 



(18) 



For a profile that is steeper than the y6-model, the exponent is 
even smaller and the A(z) dependence is even less important. 

In the case of the radial profile for entropy and pressure 
which depend on the temperature profile, we use the observa- 
tional results that entropy is approximately proportional to radius 
(e.g. Pratt 2010), which implies Tx(r) oc r~'^ 3 and the scaling for 
pressure listed in Table 1 . For Yx, which is defined as gas mass 
times global temperature and is an integral quantitiy, we take the 
global temperature to be independent of radius which then im- 
plies that Yx has the same radial dependence as gas mass and 
total mass. 

It is thus clear, that in the frame of high precision cosmology, 
we cannot just drop the variable overdensity scenario, as long as 
the recent formation approximation provides a precise picture. 
Therefore, in the next step, we critically test this approximation. 

4. Testing the density scaling with redshift by 
N-body simulations 

With the discovery that in the temperature - mass or velocity 
dispersion - mass relation the overdensity parameter can be ne- 
glected (Evrard et al. 2008), we observe a change in the literature 
in the use of scaling relations: in the recent literature a scaling 
with fixed overdensity and no overdensity evolution is prefered. 
Therefore we will in this chapter test which of these relations is 
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Fig. 4. upper panel: Testing the redshift evolution of the self- 
similar scaling of the dark matter density distribution with 
the variable overdensity scaling model. The density profiles 
are given for redshifts of z — (heavy line) and z = 
0.2521, 0.5073, 0.7695, 1.0013. Dashed lines show the un- 
sealed and solid lines the scaled profiles. Lower panel: red- 
shift evolution of the self-similar scaling of the dark matter den- 
sity distribution with the fixed overdensity scaling model. This 
model clearly works better than the model shown above. We also 
show, displaced by a factor of 0.1, the density profiles after ap- 
plying the full scaling corrections described in section 5. Apart 
from some numerical fluctuations we observe a perfect fit of the 
scaling corrections. 



better described by simulations. We are still concerned with the 
pure gravitational case and therefore apply our test to the dark 
matter distribution. Thus the most fundamental test we can per- 
form here is to study if the dark matter density distribution scales 
as Eqs. (15) and (16), which we will term the variable overden- 
sity scaling model, or just proportional to E(z) 2 , to which we will 
refer as fixed overdensity scaling model. 

For the test we use the simulations by Dolag et al. (2004) 
which are based on re-simulations of clusters taken from a large 
cosmological simulation described in Yoshida et al. (2001) and 
Jekins et al. (2001). The cosmological simulation was performed 
with 512 3 particles in a 479 lr l kpc side length box and cosmo- 
logical parameters of h — 0.7, Q„, o = 0.3, Q.a = 0.7, and a power 
spectrum normalization of <xg = 0.9. The resimulations have a 



mass resolution in the range 2 x 10 9 to 6 x 10 9 h M and a 
gravitational softening parameter of 5 h kpc. The simulation 
data set used here is for dark matter particles only to sample the 
purely gravitational evolution of the clusters. 

In Fig. 4 we show the test of the two scaling models. In the 
upper panel, which shows the variable overdensity model, we 
note that the unsealed density increases with increasing redshift, 
as expected, but when this increase is corrected by dividing by 
the factor pxA(z), we observe an overcorrection of the profiles in 
the center, where we can most sensitively test the density scaling. 
Based on all the reasoning given above, we have to conclude, 
that the central density of the clusters does not decrease as much 
with descreasing redshift as expected from the recent formation 
approximation. More explicitly, if mass accretion slowly ceases 
at low redshifts, the clusters keep a more compact shape than ex- 
pected in the simplified model. In other words, clusters observed 
at high redshift are relatively younger and more recently formed 
than cluster observed at low redshift. The recent formation ap- 
proximation is therefore not a good approximation in a ACDM 
model. 

The lower panel of Fig. 4 shows the scaling behavior of the 
fixed overdensity model. It features much better because the cor- 
rection made for the density increase with redshift is smaller and 
thus the overcorrection is less. We think that there is no more 
fundamental reason for the better match of this model than just 
a fortuitous smaller overcorrection. 

At larger radii, r > 0.15 x rsoo, the profiles rescaled with 
the variable overdensity model fit very well and actually slightly 
better than in the fixed overdensity case. 

5. Overdensity scaling based on the evolution of 
the concentration parameter 

Since also the fixed overdensity scaling model is not perfectly 
describing the evolution of the dark matter density distribution 
we now seek a more perfect description for the gravitational scal- 
ing relations and take a closer look at the predictions from N- 
body simulations. Since a lot of effort was put into an evolution- 
ary description of the universal dark matter halo mass profiles, 
the answer to our problem should be sought in this approach. 

In the literature the change of the shape of the dark matter 
halo profiles with redshift is described in the frame of the NFW 
model (see. Eq. 13) by means of the concentration parameter, 
which depends on the halo mass and the redshift of formation. 
The concentration is defined as: 



where r& is the fiducial cluster radius which can be taken as 
r A(z), ^och or r m 2oo, with ^200 being the radius for a mean over- 
density of 200 above critical density and r„,2oo refers to mean 
background density (see e.g. Duffy et al. 2008). The redshift 
evolution of this parameter has been studied by Navarro et al. 
(1997), Bullock et al. (2001), Eke et al. (2001), Dolag et al. 
(2004), Duffy et al. (2008), and Gao et al. (2008). We will use 
here the results of Dolag et al. (2004) which comes from the 
same simulations as used for our testing. The result they find is 
well approximated by 

c,„20o<*U +zT l (20) 

The same result is derived by Duffy et al. (2008, Table 1). 
Gao et al. (2008) do not give the redshift evolution explicitly, 
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but get similar results with slightly flatter evolution. The earlier 
work finds qualitatively similar results with small discrepancies 
discussed in Dolag et al. (2004). For high precision cosmology 
the current exercise should be based on simulations with higher 
statistics to be obtained in the future. But the methodic approach 
will still be the same as outlined here. 

We use the findings of Dolag et al. for the concentration pa- 
rameter to impose a second order correction to the above rela- 
tions. We apply the corrections to the relations with fixed over- 
density, since they have a simpler form and they are closer to the 
simulation results. Thus we have to transform the behavior of 
c m 200 to fixed overdensity scaling with respect to critical density, 

C= C200- 



5(Z) = C m2 Q0(z=0) (1+Z) 



-1 ^200 



r m 200 



c m200{z= o)(l+zr l n m (z) l/2 (21) 



Thus we have to correct the radius scaling by the additional 
factor. As clusters get less compact with increasing redshift, their 
radii will be larger than expected and have to be scaled it down 
accordingly by: 



r S cal = 7"200 X 



c(z = 0) 



= r a(z) 



(22) 



To observe mass conservation the change in radial scaling 
has to be compensated by a corresponding scaling of the density 
normalization. For the NFW profile the density normalization 
depends on the concentration parameter through the proportion- 
ality 



(23) 



[ln(l + c) - c/(l + c)] 



The amplitude of the density profile will be less than ex- 
pected with increasing redshift, as clusters get less compact, and 
we have to scale the density up accordingly by: 



Pscal = P X 



Pc*(z) 
PcAz = 0) 



= P/3(zT 



(24) 



The correction terms a(z) and /3(z) in Eqs. 22 and 24 are 
shown in Fig. 5 (together with the correction terms that apply 
for the recent formation approximation approach). The magni- 
tude of the correction is 35% for a and 17% for /3 at z = 1. 
The values are derived for the reference cosmology model and a 
mass independent concentration parameter c{z — 0) = 5. When 
these corrections are applied to the density profiles from differ- 
ent epochs, a perfect scaling within a few percent numerical un- 
certainty is obtained, as shown in the lowel panel of Fig. 4. 

The correction has a significant effect on the central density, 
while the effect is minor at larger radii (r > 0.2rsoo) as can be 
seen implicitly from the lower panel of Fig. 4, as the corrections 
a(z) and j3{z) have a compensating effect. In general the com- 
pensating effect for the radial profile of an observable can be 
evaluated (analogous to Eq. 9): 



Obs cor {z) = Obs(z) p(zy 



dlogObsiz) 

a(z) dlogr 



(25) 



In general we expect this correction to be small (few per- 
cent), and we will not further elaborate on this here, because 
it will be more useful once we have a better description of the 
change of the dark matter halo profiles from simulations with 
better statistics. 



n 0-8 



0.2 



redshift [1 + z] 



Fig. 5. Correction factors a(z) and (3(z) as a function of redshift 
for the fixed overdensity model (thick lines) and recent forma- 
tion approximation (rfa) model (thin lines and parameters la- 
beled with index rfa). 



6. Modifications due to hydrodynamics 

While the first part of this paper is concerned with the more the- 
oretical aspect, which model provides the best approach to de- 
scribe the scaling relations in the purely gravitational picture of 
cosmic structure growth, the second part is now exploring how 
the actually observed X-ray scaling relations of galaxy clusters 
can be described within a scenario including hydrodynamical 
effects in an empirical way. Only the most basic relations like 
those of ICM temperature or galaxy velocity dispersion with 
mass are approximately consistent with the pure gravitational 
scaling relation model as given in Table 1 . Most other relations 
involving ICM properties show deviations. The most famous of 
these deviations is the Lboi - T x relation which shows an expo- 
nent in many observational studies closer to 2.9 than to the ex- 
pected value of 2 (e.g. Edge & Stewart 1991, Ebeling et al. 1996, 
Markevitch 1998, Arnaud & Evrard 1999, Ikebe et al. 2001, Pratt 
et al. 2009). A clue to an explanation from the observational 
side comes from the fact, that the X-ray surface brightness, the 
line-of-sight integrated emission measure profile, and the density 
profile (all very closely related) can be brought to match surpris- 
ingly close at radii r > 0.15r5oo with an appropriate amplitude 
scaling (Arnaud et al. 2002, Croston et al. 2008). At smaller radii 
ICM cooling and central AGN feedback is known to modify the 
density profiles (e.g. Fabian 1994, Voit 2005). The fact that the 
shape of the ICM density profile matches so well, while the scal- 
ing of the normalization is different from the scaling in Table 1 
implies, that the ICM gas mass fraction is not constant as a func- 
tion of cluster mass. An early discussion of the change of the 
gas mass fraction with cluster mass can be found in David et al. 
(1990). 

Specifically, Arnaud et al. (2002) find a line-of-sight emis- 
sion measure scaling of the form: 



EM(r scaled ) oc T 



1.38 



instead of : EM(r sca i e d) tx T 



(26) 



where the second relation stands for the purely gravitational 
scaling. Since the integrated emission measure is proportional 
to the squared density and the line-of-sight integration path (oc 

oc r° 5 ), the amplitude scaling of the ICM density is given 
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by oc r - 44 . Croston et al. (2008) find for the matching of the 
deprojected density profiles a best fitting scaling of 

n e (r > 0.15r 500 ) * r ' 525 (27) 

for the clusters of the REXCESS sample which span a tem- 
perature range of 2 to 10 keV. Pratt et al. (2009) present di- 
rectly the gas mass fractions of the REXCESS clusters at rsoo 
and adding also the results of Vikhlinin et al. (2006), Arnaud et 
al. (2007) and Sun et al. (2009); they find a best fitting descrip- 
tion of 

/ gas( 500) <* M 02 cc r ' 3 (28) 

The variation of the temperature exponent seen in these re- 
sults is partly an effect of sample variance, but also due to the 
fact, that the ratio of the gas mass fraction of clusters of different 
mass is also a function of radius (see e.g. Pratt et al. 2010). 

Neglecting this radial dependence as a higher order effect, 
we will explore further the consequences of the mass dependent 
gas mass fraction on the scaling relations adopting a mean value 
of 0.45 for the temperature exponent of the gas fraction relation. 
Thus we find a variation of f gas with cluster mass of the form: 

M gas 045 03 

lvl tot 

The decreasing gas mass fraction with decreasing system 
mass is explained within the structure formation scenario by an 
increasing specific energy introduced into the ICM by star for- 
mation and AGN feedback (e.g. Voit 2005). In more theoretical 
approaches the reasoning for the modification of structure is de- 
rived from entropy arguments (e.g. Ponman et al. 1999, Bryan 
& Voit 2000, Voit 2005, Ostriker et al. 2005). For the present 
derivation given here we prefer the more direct observational 
approach based on the gas mass fraction. 

With this observationally implied modification to the scaling 
relations including the f gas variation with cluster mass, we can 
obtain a new set of scaling relations that are approximately con- 
sistent with the observations. These relations are given in Table 
2 (where column 2 gives the redshift evolution assuming that the 
fgas - Tx relation does not evolve with redshift). 

Compiling all the major scaling relation data from the liter- 
ature and adding new data for high redshift clusters from the 
XDCP project (Bohringer et al. 2005, Fassbender 2008) and 
from other newly detected distant clusters in Reichert et al. 
(201 1), we found first significant constraints on the reshift evo- 
lution of some major ICM scaling relations. For the mass - tem- 
perature relation a result of 

M500 = 0.29K+0.031) 7^ 62(±a08) E (z)- im(±om) (30) 

is found with mass given in units of 10 14 M Q and Tx in 
units of ke V. The relation shows a redshift evolution that is close 
to the one expected in the gravitational scenario (expectation 
E(z)~ l ). For the relation of bolometric luminosity and temper- 
ature (which is based on two independent observational param- 
eters) they find: 

Lboicsoo) = 0.079(+0.008) T 2 x 10(±0 - 24> £( z )-°- 23 (-a«) (31) 
and a mass luminosity relation of 

M 500 = 1.640*0.07) l£3£ST E(z)-°- 9 ^ (32) 



8.0 




1.0 1.2 1.4 1.7 2.1 2.5 

redshift (1+z) 

Fig. 6. Evolution parameter, E(z) 2 , as a function of redshift. The 
dashed lines show the power law functions of (1+z) with expo- 
nents of 1.26, 1.36, 1.5, 1.68 corresponding to the best approxi- 
mations given in Table 3. 



Table 3. Approximations of E(z) in powers of (1 + z) 



redshift limit power of E(z) max. deviation 



0.5 


0.63 


0.016 


0.7 


0.68 


0.027 


1.0 


0.75 


0.045 


1.5 


0.84 


0.075 



where Lboi is given in units of 10 44 erg s _1 . The last relation 
can in principle also be derived from the first two, and these 
results are consistent with each other within the error limits. 

Therefore we will focus on the implications of the luminosity 
- temperature relation. To satisfy the redshift evolution of this 
relation by means of a redshift dependent gas mass fraction, we 
have to imply a relation of the form: 

f gas oc r - 45 E(zr°- 6l5 (-°°>) (33) 

This relation can now be folded into the relations of other 
properties according to their dependence on f gas as listed in 
Table 2, to find their "empirically" predicted redshift dependence 
which is listed in column 4 of the Table. 

7. Expressions as powers of (1 +z) 

In past studies trends of cluster evolution and the evolution of 
scaling relations has most often been modeled with a redshift 
dependence of powers 1 + z in the absence of any better knowl- 
edge. For comparison with these results we studied how well the 
parameter E(z) can be expressed by powers of 1 + z. For the red- 
shift range from zero to a given upper limit we have determined 
which power of 1 + z provides the best approximation with the 
smallest maximum deviation from the function. Table 3 provides 
the results for upper redshift limits of z = 0.5, 0.7, 1.0 and 1.5. 
These approximations and the functions are shown in Fig. 6. The 
Figure immediately reveals that the approximations are not very 
good, but for limited redshift ranges (as typical for the present 
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Table 2. Hydrodynamic scaling relations: 



property 


proportionality 


scaling no evoluion off gas 


scaling with empirical evolution 


density profile 


x P fgas 


oc 


T QA5 E(z) 2 


oc r 045 ^) 1 - 385 


luminosity" 


L x <xp 2 R 3 & s 


oc 


T X A E 


oc T 2 /E-°- 23 


bolom.lum. 


L bB ,txp 2 R 3 T x ,2 f 2 as 




1 x t - 


oc Tf E- 023 


entropy, K 


oc T x p- 2 ' 3 f g 2 J 3 


oc 


T x E * /J 


oc r»- 7 £'- 74 


gas mass 


H J gas 


oc 


1 x E 


oc T x 35 E- lM5 


Yx 


OC Ty P f„n c 
A t-> j gas 


oc 


t2 95 i7-l 
X 


cc T 2 X 95 E- [M5 


pressure, P(x) 


oc Tv n f 
- 1 X J gas 


oc 


T I 45 T7l 

l x E 


oc t 1 / 5 E [M5 


surf .bright., S x(r) 


cc p 2 R f 2 

r " J gas 


oc 


T 1 / E 3 


oc A 4 E 1 77 

X 


temperature 


T x 


oc 


M 2 ' 3 E 2 ' 3 




gas mass fr. 


fgas 


oc 


M°- 3 E°- 3 


cc M°- 3 E-°- 32 


density profile 


P 


oc 


M°- 3 E(z) 23 


cc M 03 E 169 


luminosity" 


L x 


OC 


M 1.6 £ 2.6 


cc M Lb E 137 


bolom.lum. 


Lboi 


oc 


M im £ im 


cc M L93 E 2m 


entropy 


K(r) 


Cv 


M 0A7 g-0.87 


oc A/ 47 E° 3S 


gas mass 


'"gas 


OC 


M 1.3 £ 0.3 


cc M L3 E-° 32 


Y x 


T x M gas 


oc 


M 1.97 £0.97 


cc M 191 E 035 


pressure 


P(r) 


oc 


M 0.97 £2.97 


oc A/ 97 E 2 - 3S 


surf.bright. 


Sx(r) 


oc 


M 0.93 £3.93 


oc M - 93 E 2M 



Column 1 and 2 give the property and its definition, column 3 provides the scaling relation modified for hydrodynamical effects as described 
in section 6 (note that the property of column 1 is for this case assumed to be taken at the radius of the overdensity A(z)). For the redshift 
dependence noted in column 3 the f gas Tx relation of Eq. (28) is assumed to be redshift independent. For column 4 we then include the redshift 
dependence implied from the results of Reichert et al. (201 1) as given in Eq. (33). 
a for temperature independent emissivity (see Tab. 1 and text) 



surveys) the deviations are still much smaller than any typical 
observational uncertainties and therefore these simplified rela- 
tions may still be helpful at the present stage. 

Thus for the relation of bolometric luminosity and tempera- 
ture (in the fixed overdensity approximation) we find for exam- 
ple: 

Lwccrf (1+z)- 0173 (34) 
for the best fit in the redshift range 7 = 0-1. 

8. Comparison to some recent simulations 

In the Millenium gas simulations (Springel et al. 2005, Hartly 
et al. 2008) the scaling relations were studied by Stanek et al. 
(2010). For the zero redshift simulations they find slopes of the 
scaling relations with mass of 0.559(+0.002) and 0.576(+0.002) 
for the mean temperature and spectroscopic like temperature 
compared to an expectation value of 2/3. For Lt, i and Y x they 
get 1.825(±0.003) and 1.868(+0.006) compared to 2.018 and 
1.9667, respectively, as given in Table 2. The results are in fair 
agreement with the major difference that the Tx - mass relation 
is less steep in the simulations which is then also reflected in a 
slightly shallower slope of all other relations. 

For the evolution of the parameters they find that /icm, Lboi 
and Yx evolve with E(z) with an exponent of -0.44, 1 .39, and 
0.352 compared to -0.315, 1.39 and 0.33. Except for L\, a \ there 
is again a fair agreement with a slightly larger negative evolu- 
tion of the gas mass fraction in the simulations compared to the 
implication from the results by Reichert et al. (201 1). 

Using the same basic simulations Short et al. (2010) inves- 
tigated the dependence of the scaling relations on the feedback 
physics used in the simulations. As already discussed in detail 



in Reichert et al. (2011), there is a fair agreement of the obser- 
vations with the preheating models which involve an early input 
of energy and elevation of entropy of the ICM. In contrast, late 
feedback models with most of the energy input at redshifts below 
1 are clearly inconsistent with the data. 



9. Comparison to observations 

In this section we compare the model predictions to recent obser- 
vations. We use a representative set of more recently published 
results and do not aim for a complete coverage of the literature. 
In particular the earlier results either suffer from low statistics 
or they rely on lower quality observational data. A lot more re- 
sults on scaling relations and their evolution with redshift are 
expected to come in the near future with a more comprehensive 
exploitation of the XMM-Newton and Chandra archives and the 
completion of several large-survey projects, which will be used 
for a more critical test and refinement of the model. 

One of the caveats to keep in mind in the interpretation of 
the following results is, that some of the cluster surveys in the 
literature are affected by selection bias effects. These bias ef- 
fects arise e.g. from the use of flux-limited surveys, which tend 
to sample preferentially the more luminous clusters in any dis- 
tribution. We have discussed and modeled this effect for these 
types of data sets in the paper by Reichert et al. (201 1). A num- 
ber of the data sets come from the analysis of galaxy clusters 
in the data archives; also these clusters are usually studied as a 
result of discovery in flux limited surveys. Only in a few cases 
e.g. Ikebe et al. (2001) as an example for an earlier paper and 
Vikhlinin et al. (2009) have efforts been made to correct for the 
biasing effects. In most cases the bias effects on the slope of the 
relation, which is what concerns us most here, is smaller than 
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Table 4. Observationally determined slopes of the M - T relation Table 5. Observationally determined slopes of the L - T relation. 
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11 limited to systems with temperature > 3 keV 

b limited to systems with temperature > 0.6 keV 

c limited to systems with temperature > 3.5 keV 

d masses determined assuming isothermal ICM 

e core excised L x and T with r > 0.2r 50 o 

f temperatures for r = 0.2 - 0.5 x r^oo 

g core excised temperatures with r = 0.15 - 1 X r 500 

The references are from top to bottom: Finoguenov et al. 2001, Allen & 
Fabian 2001, Sanderson et al. 2003, Ettori et al. 2004, Arnaud et al. 2005, 
Maughan et al 2006, Vikhlinin et al. 2006, O'Hara et al. 2007, Hicks et al. 
2008, Zhang et al. 2008, Vikhlinin et al. 2009, Reichert et al. 2011 



L X \ is the temperature in the 0. 1 to 2.4 keV band, L X i for the 0.5 to 
2 keV band, and L] x i is the bolometric luminosity. 
a for total luminosity 

b luminosity and temperature corrected for cool core contribution 

c temperature was determined by allowing for an additional 

cool component not considered in the correlation 

d core excised L x and T with r > 0.2r 50 o 

e temperatures for r = 0.2 - 0.5 x r^oo 

f core excised L and T with r = 0.15 - 1 X r 50 o 

References not listed in Table 4 from top to bottom are: Markevitch 1998, 
Arnaud & Evrard 1999, Ikebe et al. 2002, Maughan 2007, Pratt et al. 2009 



the statistical uncertainties, and therefore we do not consider the 
selection bias effects in the following discussion . 

9.1. Mass - temperature relation 

We begin with the comparison for the mass temperature rela- 
tion. Table 4 provides a large, representative, but not complete 
list of literature results. The results by Reichert et al. (201 1) are 
based on a compilation of data from 14 data sets taken from 
the literature supplemented by recent published results on indi- 
vidual distant galaxy clusters. It gives therefore a summary or 
average of the largest and most recent observational data sam- 
ples. Most of the values for the correlation slope range from 1.5 
to 1.7 for an expected value of 1.5. The observed relation is thus 
slightly steeper but it is more close to the expectation for samples 
without low mass (low temperature) systems. For samples with 
a lower temperature limit above 3 keV the slopes are shallower 
and closer to the self-similar scaling. 

9.2. X-ray luminosity - temperature relation 

The luminosity temperature relation is the relation with the two 
observables derived almost independently. The luminosity is ob- 

5 In most cases, where the cluster sample has been compiled from 
data archives and not from very well defined surveys with published 
selection criteria and sensitivity functions, a rigorous reconstruction of 
the selection effects is not possible anyway 



tained from imaging data with tiny corrections from spectral 
data, while the temperatures originate from the interpretation of 
X-ray spectra. In Table 5 we summarize the observational results 
for the luminosity - temperature relation for both, luminosities 
derived for certain energy bands and bolometric luminosities, 
Lboi- 

The band limited luminosity scaling relations have observed 
slopes in the range 2 - 2.5 with an expected value of about 2.4, 
while most of the Lboi - T relation slopes show observed values 
of 2.6 - 3.7 with an expected value of about 2.9. The results 
are thus in good agreement with the predictions of the modified 
scaling relations within the observational uncertainties. 

9.3. Luminosity - mass relation 

The X-ray luminosity - mass relation is one of the most impor- 
tant relations for cosmological modeling of X-ray cluster sur- 
veys. Published results are listed in Table 6. Again we find good 
agreement of the observed slopes of 1.4 - 1.7 (for Lb an d) and 1.6 
- 2 (for Lboi) with the predictions for the slope of 1.6 and 1.93, 
respectively. 

9.4. Entropy - temperature relation 

Table 7 lists the results for the entropy temperature relation from 
the observational analysis of Pratt et al. (2010). We note that 
the slope of the relation depends on the scaled radius at which 
the measurement is taken. In contrast to the observational pa- 
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Table 6. Observationally determined slopes of the L - M relation. 



relation 


slope 


comments 


reference 


Lxi 


■'"200 


1.61 ±0.09 


106 cluster, ROSAT/ASCA 


Reiprich02 


Lxi 


- M a) 

200 


1.46 ±0.11 


63 cluster, ROSAT/ASCA 


Reiprich02 


Lboi 


- M"' 
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lUo cluster, KUoAl /A^CA 


Reiprich02 


i 
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- M 20 o 
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1 1 Hunters 7 — f) fS — 1 


lVTan fftinnOfS 


I , , 
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~ M 500 
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Maughan07 
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1.45 ± 0.07 




Maughan07 


Lbol 


- M soo 


1.03 ±0.28 


13 RCS Clusters, z = 0.6 - 1.1 


Hicks08 


Lbol 


-M 500 


2.33 ± 0.70 


37 clusters, XMM - Newton 


Zhang08 
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- M 500 


1.61 ±0.14 


17 cluster, CHANDRA 


Vikhlinin09 
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1Y1 Y3O0 
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1Y1 Y3O0 
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- M d) 
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1.90 ±0.11 
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ArnaudlO 
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-M d) 
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1.76 ±0.13 




ArnaudlO 


Lboi 


-M 500 


1.51 ±0.09 


14 literature samples 


Reichertl 1 



L X i is the temperature in the 0. 1 to 2.4 keV band, L X i for the 0.5 to 

2 keV band, and Lboi is the bolometric luminosity. 
a Quoted is the BCES orthogonal fit. 
b Core excised luminosity, r = 0.15 - 1 x r^oo 
c My is the mass estimated from the Y x - M relation. 
d Corrected for Malmquist bias 

References not listed in Tables 4 and 5 from are: Reiprich & Bohringer 

2002, Arnaudet al. 2010 



Table 7. Observationally determined slopes of the Entropy - T 
and entropy - mass relation 



relation 


slope 




comments 


reference 


A'soo - T 


0.92 ± 


0.24 


31 clusters XMM 


-Newton PrattlO 


K iooo _ T 


0.83 ± 


0.06 




Pratt 10 


K 2500 - T 


0.76 ± 


0.06 




PrattlO 


^1000 ~~ T 


0.83 ± 


0.06 




PrattlO 


K$oo - M500 


0.62 ± 


0.17 




PrattlO 


^2500 _ -^2500 


0.42 ± 


0.05 




PrattlO 



11 limited to systems with mass > 3 x 10 13 M G 
b limited to systems with temperature > 0.6 keV 
The references is Pratt et al. 2010 



rameters above which are integrated values like e.g. luminos- 
ity, this is the local entropy value at a given radius. The change 
in slope originates from the fact that the gas density reduction 
in low mass systems is not completely self-similar with radius, 
somewhat contrary to what is suggested e.g. by the very tight 
scaling of the density profiles in Fig. 2 of Croston et al. (2008). 

The gravitational scaling model would predict the relation 
K oc T, while our modified scaling relation predicts K oc T 01 . 
The latter prediction is approximately met at intermediate radii. 
For larger radii the results seem to approach the gravitational 
scaling relations. 



Table 8. Observationally determined slopes of the M gas - T rela- 
tion 



relation 




slope 




comments 




reference 


M^,200 - 


T 


2.22 ± 


0.31 


11 clusters, z = 0.6 


- 1.0 


Maughan06 


Mg,2500 - 


- T 


1.80 ± 


0.49 






Maughan06 


M giQ0 ~ 
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2.12 ± 


0.12 


31 lowzcl., XMM 


- Newton 


Croston08 


^,500 - 


-pa) 


1.86 ± 


0.19 


37 clusters, XMM - 


- Newton 


Zhang08 



il Temperature for r = 0.2 - 0.5 x 7-500 

The reference not listed in previous Tables is Croston et al. 2008 



9.5. Gas mass - temperature relation 

Some results for the gas mass - temperature relation are listed 
in Table 8. The observed range of values for the slope of 1.88 - 
2.22 corresponds well to the predicted value of 1.95. 

9.6. Other relations 

Observational results for some other interesting relations are pre- 
sented in Table 9. The Yx parameter is a quantity motivated 
by cluster observations in the Sunyaev-Zeldovich effect (SZE). 
Since the SZE signal is proportional to the intracluster plasma 
temperature and the total number of electrons, the Yx-parameter 
is defined as Yx = T X M gas . This parameter has been promoted 
as a very good mass proxy for galaxy clusters (Motl et al. 2005, 
Kravtsov et al. 2006). 

Values for the slope of the L\,„i - Yx relation of 0.65 - 1.1 
and for the Lb an d - Yx relation of 0.8 - 1.14 agree well with the 
predictions (combining the Lx-T and Yx— T relations) of values 
of 0.98 and 0.81, respectively. The M500 - Yx relation with slope 
values of 0.57 - 0.62 is also not far off from the precdited value 
of 0.51. 

Finally the M tot —M gas relation with a predicted slope of 0.76 
is shallower than the listed result of 0.9 by Zhang et al. (2008) 
and also somewhat shallower than the slope inferred from Fig. 9 
in Vikhlinin et al. (2009) which roughly corresponds to a slope 
of 0.83 in the mass range 10 14 to 10 15 M . But our prediction is 
very well supported by (and partly originates from) the results 
of Pratt et al. (2009) based on a large and representative sample 
of galaxy clusters with high quality XMM-Newton data. 

9.7. Evolution with redshift 

Vikhlinin et al. (2009) quote a result for the evolution of the 
Lxi - M relation of 

L X2 oc M 161±014 £ (z )l-85*0.42 (35) 

which has to be compared to Eq. (32) based on the result 
of Reichert et al. (2011). Inverting relation (32) and converting 
from L bo i to L x2 wefmd:L z2 <* L bol T- i)5 oc T- Q - 5 M L92 E(z) i n . 
Substituting Eq. (7) without the A-factor we get: 

Lxi <* M U9 (S) £( z ) L40 (-a7) (36) 

showing agreement between the two results within the un- 
certainties. 

More effort has been put on the study of the evolution of the 
luminosity - temperature relation. Here we find a different pic- 
ture, such that most of previous works finds a positive evolution. 
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Table 9. Observationally determined slopes for various relations 



relation 


slope 




comments 


reference 
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y x 


1.10 ± 
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Hicks08 
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Zhang08 
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0.906 + 0.08 




Zhang08 


M 500 ■ 
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Zhang08 
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0.57 ± 
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17 clusters, CHANDRA 
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0.84 ± 
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1.07 ± 
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a Core excised luminosity, r 
b core excised L and T with r = 
c Corrected for Malmquist bias 



0.15 - 
= 0.15 



1 x r 500 
- 1 x r 5 oo 



In the past literature, the evolution has most often been parame- 
terized in the form: 



T p (1 +zf 



(37) 



Taking this functional form and the correspondence of 
E(z) ~(1 + z) 75 from Table 3 for the redshift range z = — 1 
we find for the results of Reichert et al. (201 1) a proportionality 
of the evolution of the Lx - T relation of oc (1 + 2) _017 (-o:«,). 
Earlier works find a — 1.5 + 0.3 in Vikhlinin et al. (2002) and 
Lumb et al. (2004), and a = 1.8 ± 0.3 in Kotov & Vikhlinin 
(2005). Maughan et al. (2006) find a value of or = 1.3 ± 0.2 if 
they combine their WARPS cluster sample with that of Vikhlinin 
et al. (2002) and a value of or = 0.8 ± 0.4 for the WARPS sam- 
ple alone. Ettori et al. (2004) find the least positive evolution 
with values for the exponent of (1+z) in the range 0.04 to 0.98 
when using Markevitch (1998) and Arnaud & Evrard (1999) as 
local reference, but partly negative values in the range -0.48 to 
+0.54 when combining their results with those of Novicki et al. 
(2002). Common to all these studies is their use of the data of 
Markevitch (1998) and Arnaud & Evrard (1999) for the local 
(z ~ 0) reference. The normalization of the Lx - T relation in 
the latter two works is lower by a factor of very roughly 1 .5 than 
that of Pratt et al. (2009) which has been used in Reichert et al. 
as the most important local reference. Therefore, the most ob- 
vious reason for the fact that these earlier works find a positive 
evolution in contrast to the slightly negative evolution found by 
Reichert et al. (201 1) is the different local reference in addition 
to the small sample sizes and selection bias effects. 

Pacaud et al. (2007) also quote an Lx - T evolution best fit 
by a — 1.5 + 0.4, but the fit is poor and the data are different 
from this approximation at z > 0.7 as can be seen in their Fig. 
4. Maughan et al. (2011) in a recent paper also find a positive 
evolution with substantial deviations from the overall trend at 
intermediate redshifts. Different from these other results O'Hara 
et al. (2007) find a negative evolution with a value of a — -0.25+ 
0.56. 

In summary, it is clear that the sparsity and inhomogeneity 
of the data used in the past to study the redshift evolution of 
the scaling relations lead to inconclusive results and we are just 



beginning to see some trends now. The uncertainties of the pa- 
rameterized evolution of the relations are still very large and e.g. 
in the case of the Lx - T relation the results are still consistent 
with no evolution of the Lx - T relation. 



10. Discussion and Conclusion 

Studying in detail the evolution of the dark matter mass den- 
sity profiles of simulated galaxy clusters, we have shown that 
the older model of self-similar scaling relations based on the 
recent formation approximation which uses a scenario where 
the fiducial overdensity radius, r^ z ), is taken to be redshift de- 
pendent, is not accurate and a scaling with a fixed overdensity 
provides a better and currently sufficiently precise description 
of self-similar evolution. For a precise analysis of future cos- 
mological surveys of the galaxy cluster population, we should 
improve this description further. We plan to do this with larger 
N-body/hydrodynamical simulations which are performed at 
present and therefore give no detailed recipes for the scaling cor- 
rections in this paper. The corrections as shown in Fig. 5 depend 
on the cosmological model used and therefore these corrections 
will have to be worked separately for each model case studied. 

Studying the different types of scaling relations involving pa- 
rameters derived from X-rays, we can distinguish two types of 
scaling behavior: the mass - temperature relation is mostly de- 
pendent on the scaling of the dark matter potentials and is there- 
fore very close to the gravitational scaling prediction. Most other 
relations involving gas density or gas mass are affected by the 
non-constancy of the gas mass to total mass ratio. With the intro- 
duction of modified scaling relations to take this hydrodynamical 
effects into account, we can describe the currently available data 
sets within the given uncertainties. 

Looking at the evolution of the scaling relations with red- 
shift, we find an analogous situation: the M - T relation corre- 
sponds within the current uncertainties to the prediction of the 
gravitational self-similar scenario. All other relations involving 
parameters which depend on the gas density show deviations, 
which implies that the gas mass fraction is not constant for given 
cluster mass with redshift. This is just the consequence of the 
following effects. At higher redshifts clusters of given mass are 
more compact and the ICM has to be squeezed into a deeper 
and narrower potential. Since in preheating models, which seem 
to explain the data best, the gas starts out with an elevated en- 
tropy before cluster formation is complete, the gas is less tightly 
squeezed into the earlier, narrower potentials than into the later 
wider potentials. 

Our modified scaling relation model does not describe all ob- 
servational effects. As shown in Pratt et al. (2010), the entropy 
scaling depends on the radius at which the entropy is measured. 
These results imply that the ICM depletion is larger in the center 
of groups and clusters than in the outer parts. More data are re- 
quired that extend out to large cluster radii (to rsoo and beyond) 
to substantiate this result. 

The redshift evolution of the L — T scaling relation, which is 
predicted to be oc E{z) in the simple self-similar scenario, is now 
found in recent studies to be much less positive or even negative 
(e.g. O'Hara et al. 2007, Reichert et al. 2011). This has the im- 
portant consequence that one will find less high redshift galaxy 
clusters in future X-ray and SZE surveys, than predicted based 
on the simple scaling models (Reichert et al. 201 1). More impor- 
tantly, a precise measurement of this evolution effect is crucial 
for using the future X-ray survey data on galaxy clusters for the 
test of cosmological models. 
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As shown by the comparison of the evolution of the scaling 
relations compiled by Reichert et al. (201 1) and the simulations 
by Short et al. (2010), the study of the evolution of the scaling 
relations also provides important insight into the astrophysics 
of the ICM. Currently the observational data strongly favour a 
model with ealry preheating of the ICM. 

In the coming years both the observational data as well as the 
simulation results will experience further strong improvements. 
Therefore we see the importance of this paper more in elucidat- 
ing the way how the scaling relations should be analysed and 
applied, rather than already providing the best parameterization 
of the results. 
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